Disorder-induced phases in the 5 = 1 antiferromagnetic Heisenberg chain 
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We use extensive density matrix renormalization group (DMRG) calculations to explore the phase 
diagram of the random 5 = 1 antiferromagnetic Heisenberg chain with a power-law distribution of 
the exchange couplings. We use open chains and monitor the lowest gaps, the end-to-end correlation 
function and the string order parameter. For this distribution, at weak disorder the system is in 
the gapless Haldane phase with a disorder dependent dynamical exponent, z, and z — 1 signals 
the border between the nonsingular and singular regions of the local susceptibility. For strong 
enough disorder, which approximately corresponds to a uniform distribution, a transition into the 
random singlet phase is detected, at which the string order parameter as well as the average end-to- 
end correlation function are vanishing and at the same time the dynamical exponent is divergent. 
Singularities of physical quantities are found to be somewhat different in the random singlet phase 
and in the critical point. 

PACS numbers: 64.60.Ak, 05.50.+q, 68.35.Rh 



I. INTRODUCTION 

The S = 1 spin antiferromagnetic Heisenberg chain has 
received much attention, both experimentally and theo- 
retically, since Haldane 1 conjectured that its low-energy 
properties are qualitatively different from that of the ex- 
actly solved 5 — 1/2 model. The 5 = 1 chain (together 
with all other integer spin chains) has a finite gap in the 
excitation spectrum and hidden topological order in the 
ground state, which is characterized by the string corre- 
lation function^ On the other hand, the bulk spin-spin 
correlations of the model are short ranged, having a fi- 
nite correlation length, £. In an open chain of length L, 
there are spin 5 — 1/2 degrees of freedom at each edge 
and the end-to-end correlations approach a finite value 
in an exponential fashion, having the same characteristic 
length scale, £, as bulk correlations. 3 

Quenched disorder, which is realized by random cou- 
plings, also has different effects for 5 = 1/2 and S = 1. 
In the former case any amount of disorder is enough to 
drive the system into a new type of fixed point £ whereas 
for the S = 1 chain, weak disorder is irrelevant and the 
properties of the weakly random chain are the same as 
that of the pure one£ For stronger disorder, however, 
the low-energy properties of the system are changed and 
detailed analytical and numerical investigations were de- 
voted to clarify the properties of the new random fixed 
points. 

The analytical studies of the random chain are made 
by variants of the strong disorder renormalization group 
(RG) method, which has been introduced for the 5=1/2 
chain by Ma, Dasgupta, and Hu& and has been ana- 
lyzed in great details by Fisheri* This strong disorder 
RG method has been used afterwards for a large variety 
of random quantum and classical systems, (for a review, 



see Ref. 0). For the 5 = 1 chain, extensions of the orig- 
inal Ma-Dasgupta rules are necessarj^S to describe the 
disorder induced phases in the system, which include a 
gapless Haldane (GH) phase, for intermediate disorder, 
and a random singlet (RS) phase, for stronger disorder. 

Numerical studies of the random 5=1 chain have 
been made by exact diagonalizationiS of the density 
matrix renormalization 11 (DMRG), by quantum Monte 
Carlo (QMC) methodsjiSii 3 . and by numerical implemen- 
tations of the strong disorder RG method^ Despite con- 
siderable numerical effort, several aspects of the low- 
energy properties of the random 5=1 chain are still 
unclear and some numerical results are conflicting. In 
the numerical calculations mainly boxlike distribution of 
disorder is considered, which, as noted in Ref. lla, repre- 
sent only a limited strength of randomness. In numerical 
RG studies both the GH and the RS phases are identi- 
fied; however, the transition point between these phases 
is rather approximate. In DMRG calculation (see also 
Ref. ITo|) Hidaii has identified only the GH phase, and 
conjectured that the RS phase is not accessible for any fi- 
nite strength of disorder. In a comment to Hida's work^i 
Yang and Hyman 15 have predicted the appearance of the 
RS phase for some type of power-law distribution of the 
disorder. Another numerical work by QMC simulations 12 
has shown the existence of the RS phase even for the 
boxlike distribution and these results are confirmed by 
independent QMC simulationsii 3 - In the QMC calcula- 
tions, some properties of the RS phase are verified (cf. 
scaling relation between length and time, decay of the 
string correlation function), but results about the spin- 
spin correlation function are different from the RG pre- 
dictions. At the critical point no numerical estimates are 
available to check analytical RG predictions. We note on 
recent studies of Griffiths effects^ in the system with en- 
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forced dimerizatio n-7^§ and related work on the random 
5 = 3/2 and higher spin chainsi2i22i2i. 

In this paper, our aim is to study the low-energy 
properties of the random S — 1 chain by the DMRG 
method^ The features of our study are the following: 
(i) We consider a more general (power-law) distribution 
of disorder, which allows us to enter more deeply into the 
RS phase, thus to obtain convincing evidence of its exis- 
tence, (ii) We calculate a different physical quantity, the 
end-to-end correlation function, which carries important 
information about the phases of the system. The aver- 
age end-to-end correlation function has a finite limiting 
value in the GH phase and vanishes in the RS phase. Fur- 
thermore, in the GH phase from the low-value tail of its 
distribution, independent estimates about the dynamical 
exponent are obtained, (iii) We try to perform a com- 
parative analysis between the properties of the system at 
the critical point and in the RS phase and to check the 
available RG predictions. 

The structure of our paper is the following. The model, 
the basic ingredients of the strong disorder RG methods, 
and the conjectured phases are given in Sec. [H] Re- 
sults of our DMRG studies are presented in Sec. IIHI and 
discussed in Sec. IIVI 



II. 



THE MODEL AND THE STRONG 
DISORDER RG RESULTS 

A. Model 



We consider the spin 5 = 1 random antiferromagnetic 
Heisenberg chain with the Hamiltonian 



H — JiSi ■ iSj+i, 



(1) 



where the J$ > are independent and identically dis- 
tributed random variables. Here, we use the following 
power-law distribution 



ps(J) = 6- 1 J- i+ V s 



for < J < 1, 



(2) 



where 5 2 = var[ln J] measures the strength of disorder. 
In previous numerical work, a boxlike distribution was 
used, 



Pw(J) 



1/W for 1 - W/2 < J < 1 + W/2 
otherwise, 



(3) 



in which the strength of disorder grows with W. Note 
that the possible maximal value, W = 2, corresponds to 
the uniform distribution, which can be obtained from Eq. 
(J2J with 5 = 1 and having a prefactor, 1/2, and a range 
< J < 2. Consequently, the power-law distribution for 
6 > 1 represents a disorder, which is stronger than any 
boxlikc disorder. 

The low-energy behavior of the system of size, L, is 
encoded in the distribution of the lowest gap, A, denoted 



by Pl(A). We note that for an open chain the first gap 
corresponds to the localized edge states; therefore, one 
should study the second (not localized) gap. The average 
spin-spin correlation function is denoted by 



(4) 



where [• • -] av stands for averaging over quenched disor- 
der. For bulk correlations with \j — i\ -C i,j = O(L), we 
have C(i,j) = Cb(\j — whereas for end-to-end corre- 
lations, C(l, L) = Ci(L). The string correlation function 
of the model is defined b^ 



O z {r) = -(Sf exp [in (S? +1 + Sf +2 + ■■■ + 5f +r _0] S? +r ) 

(5) 

and its large r limiting value is the string order parame- 
ter. For several quantities it turned out useful to consider 
the average of its inverse. More precisely, for a physical 
observable, /, we denote by / lv the following quantity: 



r = 



(6) 



[f-\ v ' 

what we shall call as inverse average. 

B. Weak disorder limit— Haldane phase 

In absence of randomness (Jj = J) the spectrum 
is gapped^ and bulk spin-spin correlations are short 
ranged, C (r) ~ exp(— r/£) with £ = 6.03\ On the 
contrary, end-to-end spin-spin correlations and the string 
correlation function have a finite limiting value. For weak 
disorder, when the distribution of J is sufficiently nar- 
row, the Haldane gap is robust and the system stays in 
the Haldane phase. 5 The border of the Haldane phase 
can be estimated by noting that the Haldane gap is ro- 
bust against enforced dimerizationjSi when even and odd 
couplings are different, so that 



Ji = J(l + £)(-!)*) exp(*Ci) 



(7) 



where Qi are random numbers of mean zero and variance 
unity. The pure system [5 = 0) for D < 0.25 stays in the 
Haldane phase— and at the phase transition point the 
coupling at an odd bond, J Q , and that at an even bond, 
J e , are related as J = 0.6J e . We expect that in the 
presence of disorder the Haldane gap stays finite, if the 
maximum (J m ax) and the minimum (J m in) values of the 
couplings satisfy J mini Jmax > 0.6. From this argument 
we obtain for the border of the Haldane phase for the 
box distribution Wg ~ 0.5. On the other hand, for the 
power-law distribution in Eq. (J2J Jmin = 0; therefore, for 
any 8 > the Haldane phase is expected to be destroyed. 



C. Strong disorder limit— RG approach 

For strong disorder the low-energy properties of the 
system are explored by variants of the strong disorder 
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RG approach. In the standard Ma-Dasgupta-type RG 
approach, the couplings of the random antifcrromagnetic 
Heisenberg chain are put in descending order and the 
largest coupling defines the energy scale, Q, in the sys- 
tem. During renormalization the pair of spins with the 
largest coupling, say Ji = CI, are replaced by a singlet 
and decimated out. At the same time a new coupling 
is generated between the spins at the two sides of the 
singlet, which is given in a perturbation calculation as 



~ 4 J»_! J i+1 
3 J t 



(8) 



As noticed by Boechat, Saguia, and Continentino^ for 
weak disorder some of the generated new couplings can 
be larger than the energy scale, O. Therefore, the stan- 
dard strong disorder RG approach works only for strong 
enough disorder and describes only the RS phase of the 
system. 

To cure this problem, different types of RG approaches 
are proposed. Monthus et alM> suggested to replace the 
pair of spin S = 1 connected by the strongest bond by 
a pair of S = 1/2. In this case the renormalized sys- 
tem consists of a set of spin S = 1 and S — 1/2 degrees 
of freedom, having both antiferromagnetic and ferromag- 
netic couplings. The renormalized couplings, which are 
calculated perturbatively, are all smaller than fl. This 
RG approach, during which no spin larger than S = 1 is 
generated, can be used to describe both the gapless Hal- 
dane and the RS phases and provides precise numerical 
estimates about the critical exponents. 

In another modified RG approach, Saguia et al*^ use 
the standard perturbative approach in Eq. © , provided 
max(Ji-i, Ji+i) < 3S1/4. Otherwise the triplet of spins 
with couplings, max(Ji-i, Ji+\), and £1 is replaced by a 
single spin. Also in this method the variation of the en- 
ergy scale is monotonic: the generated two new couplings 
are both smaller than VI. 

Recently, a variant of the strong disorder RG method 
was proposed by one of usj2& in which the pair of spin 
with the strongest coupling is decimated out, but-and 
this is a feature of our current method-the new coupling 
between the remaining spins is calculated nonperturba- 
tively. The four spins with couplings J<_i, Jj, and Jj+i 
are replaced by two spins and during decimation the low- 
est gap in the two systems remains the same. It is easy 
to see that the rule we use is somewhat similar to the 
approach by Saguia et alw^ However, this method has 
no discontinuity in the approximation, which could be 
important in the vicinity of the critical point, where a 
crossover takes place between the different decimation 
regimes. 



D. Disorder-induced phases 

Based on a modified strong disorder RG approach&2ii4 
and different numerical calculationsjii*i2ii^ the following 



scenario of the phase transition in the model is conjec- 
tured with increasing strength of disorder. 



1. Gapless Haldane phase 

For sufficiently strong disorder (8 > 8q or W > Wg), 
the gap in the Haldane phase is closed and one arrives 
to the gapless Haldane phase. As we have argued in Sec. 
Ill Bl Sq — and Wg ~ 0.5. The GH phase is a quantum 
Griffiths phased in which the correlation length, £(<5), 
is finite, whereas the typical time scale, t r ~ A -1 , is 
divergent. Relation between the size of the system, L, 
and the smallest gap is given by 

A ~ L~ z , (9) 

where z is the disorder induced dynamical exponent. The 
distribution of the lowest gap is given by 



P L (A)dA = L- Z P 



A 

17- 



dA, 



(10) 



and P{x) ~ x~ 1+1 l z for small x, so that from the low- 
energy tail z can be calculated. Similarly the distribution 
of the end-to-end correlation function has a vanishing 
tail, which behaves asSI P(C\) ~ C 1 _1+1 ^, which gives 
an independent way to calculate the dynamical exponent. 
Some thermodynamical quantities such as the local sus- 
ceptibility, x, and the specific heat, c v , are singular at 
low temperature^ 



X(T) ~ T- 1+1/z , c v {T) ~ T 1/z 



(11) 



The limit of divergence of x(T) is signaled by z = 1, 
and the corresponding disorder is denoted by 8\ (Wi). 
The separation of the two parts of the GH phase with 
z < 1 and z > 1 can be located by considering the inverse 
average of the gap, A lv , and the inverse average of the 
end-to-end correlation function, C{ v . In the nonsingular 
region, z < 1, both A lv and C{ v are finite, whereas in the 
singular region, z > 1, both are vanishing. 

To see this, we consider the inverse average of the gap 



A -2+l/ ZdA 



-1 + 1/2 



v -1+1/2 

*max 



(12) 

which indeed tends to zero, if z > 1 and A min — * 0. On 
the other hand for the vanishing of the average gap one 
needs A max — > 0. One can use a similar reasoning for 
the end-to-end correlation function, for which the upper 
limit of the distribution, C™ ax > 0, thus [Ci]; v > 0, in 
the whole region, 8 < Si. 

In a static sense, the gapless Haldane phase is noncrit- 
ical: the average end-to-end correlation function, as well 
as the string order parameter, is finite in the complete 
GH phase. 
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2. Critical point 

Increasing the strength of disorder over a critical value 
(5c or Wc), the system arrives at the random singlet 
phase. As the critical strength of disorder is approached, 
the correlation length diverges: £ ~ (6c — 5)~ u , with 
v = (1+ \/l3)/2 and the string order parameter vanishes 
aMO z (8)~ (Sc-S) 2 ?, with/3= [2(3-V5)/(Vl3-l)]. 
At the critical point the string order-parameter decays al- 
gebraically, O z (r) ~ r - ' 7 , with 77 s * = 2j3/v, The end-to- 
end correlation function goes to zero algebraically, too, 
Ci(L) ~ L - * 71 , similarly to the bulk spin-spin correla- 
tion function, Cf,(r) ~ r _?) . Here, however, there are 
no theoretical conjectures about the exponents rji and r\. 
The relation at the critical point between the correlation 
length and the relaxation time is strongly anisotropic, 



In t r 



(13) 



with ip = 1/3; thus, the dynamical exponent, z, is for- 
mally infinity. This type of infinite disorder scaling is 
seen in the distribution of the gaps, which is given by 



Pi(A)dA = L~^P 



In A 



d In A 



(14) 



In the space of variables, dimerization (D) and disorder 
(6), the critical point of the system represents a mul- 
ticritical point in which three Griffiths phases with dif- 
ferent symmetry meet^i The corresponding exponents 
follow by permutation symmetry and the calculation can 
be generalized for higher values of SH 



3. Random singlet phase 

For a disorder 5 > 5c (W > Wc), the low-energy 
behavior of the system is controlled by an infinite dis- 
order fixed point and the system is in the RS phase. 
The RS phase is a critical phase, both £ and t r are di- 
vergent, and its properties are assumed to be identical 
to the RS phase of the random S = 1/2 chain. This 
latter system is studied in great detail by Fisher— with 
the asymptotically exact strong disorder RG method and 
these results have been confronted with detailed numer- 
ical investigations. 28 29,30 Here we repeat that in the RS 
phase there is infinite disorder scaling, so that relations in 
Eqs. I|13|) and l|14|) are valid with an exponent, ip = 1/2. 
The RS phase is instable against enforced dimerization, 
as given in Eq. Q , and the correlation length behaves as 
£(D) ~ D~ u , with v = 2. In the RS phase the bulk and 
end-to-end correlation functions decay algebraically. In 
Table [I] we collected the conjectured values of the critical 
exponents both in the random singlet phase and at the 
critical point and compared these values with the esti- 
mates obtained in this paper. 



TABLE I: Theoretical predictions for the critical exponents in 
the random singlet (RS) phase and at the critical point (CP). 
Values obtained in this paper are given in square brackets. 





n st 


P 


V 


4, 


V 


m 


RS 


0.382(0.41(4)] 




2 


1/2 (0.45(5)] 


2 


1 (0.86(6)] 


CP 


0.509(0.39(3)] 


0.586 


2.30 


1/3 (0.35(5)] 




- (0.69(5)] 



TABLE II: Numerical estimates for the borders of the differ- 
ent phases of the random S = 1 chain with boxlike disorder, 
see Eq. Q. The different ph ases are defined by W < Wg, 
Haldane phase; Wg < W < W\ , GH phase with nondivergent 
local susceptibility; Wi < W < Wc, GH phase with divergent 
local susceptibility; and W > Wc, RS phase. 
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Wi 


W c 


RGf9] 






1.48 


RG[14] 




0.76 


2. 


DMRG Hi] 




1.8 


no 


QMG12] 


1.37 


1.7 


1.8 


QMCfl3] 






< 2 



E. Summary of the existing numerical results 

In previous numerical studies the box distribution in 
Eq. has been used. In Table ITT1 we present the esti- 
mates of the borders of the different phases obtained by 
different numerical methods, such as by numerical im- 
plementation of the strong disorder RG, by DMRG, and 
by QMC. We note that for the power-law disorder in Eq. 

the critical disorder is estimated^ as 5c ~ 1-5. Using 
variants of the strong disorder RG method(2iI4 the calcu- 
lated critical exponents in the RS phase-within numerical 
precision-correspond to the predicted, analytical values. 
To reach the asymptotic region, however, one often needs 
to treat very long chains of length L = 10 4 — 10 6 , see Ref. 
9. The DMRG and QMC investigations have led to dif- 
ferent conclusions for the strongest box disorder, with 
W = 2. No RS phase is found by DMRG,ii whereas 
by QMC infinite disorder scaling is detect edi 1 ^* 1 ^ The 
average string correlation function was shown to decay 
algebraically withA 3 i] s t = 0.378(6), close to the theo- 
retical result in Table HJ On the other hand, the aver- 
age spin-spin correlation function was found to have an 
exponent^ r\ = 1, which greatly differs from the theo- 
retical value of 77 = 2. 



III. NUMERICAL INVESTIGATIONS 

A. The DMRG method 

Most of our numerical results are based on DMRG cal- 
culations. In this case we used open chains up to length 
L = 64, for weak disorder and up to L = 32 for strong 
disorder and calculated the lowest two gaps, the end- 
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to-end correlation function and the string order param- 
eter. This latter quantity is calculated for open chains 
between points i = L/A and j = 3L/4. Note that for an 
open chain the first gap is related to the surface degrees 
of freedom and goes to zero exponentially with L. The 
characteristic bulk excitations are given by the second 
gap and we studied this quantity. In the numerical cal- 
culation we have retained up to m = 180 states in the 
DMRG and checked that convergence of the numerical 
results is reached. We used the power-law distribution of 
disorder in Eq. |J5| in the canonical ensemble, i.e., there 
was no constraint to the value of the sum of the odd and 
even couplings. In this way there is a nonzero residual 
dimerization, which could be the source of some error 
for small systems. However, using the microcanonical 
ensemble, in which the sum of the odd couplings is the 
same as that of the even couplings, could leave to dif- 
ferent finite-size exponents for the end-to-end correlation 
function, which is known for the random transverse-field 
Ising chaini 31 i 32 i 33 We have calculated typically 10 000 
independent disorder realizations in each case. 
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FIG. 1: The inverse average gap, A lv , as a function of the 
inverse size of the system for different strengths of disorder, 
S — 0.1,0.2,0.3,0.4, and 0.5 from the top to the bottom, 
respectively. A lv seems to vanish around 5 ~ 0.5. Inset: the 
inverse average end-to-end correlation function as a function 
of inverse size, with the same values of disorder, as in the 
main panel. Note that for weak disorder the size dependence 
of Cf is nonmonotonic, which is due to a finite correlation 
length in the system. 



B. Gapless Haldane phase 

1. Nonsingular region: z < 1 

We have calculated the distribution of the (second) 
gap and determined its inverse average, A IV , which is 
presented in Fig. El as a function of the inverse size, 
.L -1 , for different values of 8. The limiting value as L — > 
oo is monotonously decreasing with 8 and A IV seems to 
approach zero at a limiting disorder 8\ ~ 0.45 — 0.5. 
Similar conclusion is obtained from the behavior of the 
inverse average of the end-to-end correlation function, 
C{ V (L), which is shown in the inset of Fig. 

Note, however, that the average end-to-end correlation 
function, as shown in Fig. [21 is finite at S±. The extrap- 
olated values of A IV and C{ V (L) are shown in Fig. El 
Close to Si, both curves are compatible with an approx- 
imately linear variation with S± — 8. At the boundary 
point, Si = S, the size dependences of A lv and Ci{L) 
are shown in the inset of Fig. Both are linear in L _1 , 
in accordance with the criterion that at Si the disorder 
induced dynamical exponent is z(8i) = 1. 



2. Singular region: z > 1 




FIG. 2: Extrapolated values of the inverse average gap and 
the inverse average end-to-end correlation function as a func- 
tion of the strength of disorder. At S = we obtain estimates 
for the nonrandom model, A = 0.4105(3) and Ci = 0.283(1). 
Inset: Size dependence of the inverse average gap and the in- 
verse average end-to-end correlation function in a log-log plot 
at the boundary of the gapless Haldane phase. Both lines 
have an approximate slope, z = 1, denoted by broken lines. 
The typical value of the error is indicated; otherwise, the error 
is smaller than the size of the symbol. 



We have calculated the average string order parame- 
ter and the average end-to-end correlation function for 
different sizes L. Also we have determined the disorder 
induced dynamical exponent, z, which is deduced from 
the low-energy tail of the gap distribution [see Eq. itTUI) ]. 
The extrapolated values of O z and Ci , as well as 1/z are 
plotted in Fig. 13 for different strengths of disorder. All 
these three quantities tend to zero around the same limit- 
ing value of disorder and the border of the Griffiths phase, 



i.e., the location of the critical point can be determined 
as 8 C — 1.0(1). In the extrapolation procedure we have 
made use of the finite-size dependence of O z (L) ~ L~ 7,st 
and Ci(L) ~ L~ m at the critical point, which is shown 
in the inset of Fig. For weaker disorder, S < S c , O z 
tends to a finite limiting value, which is illustrated in Fig. 
01 using a scale L _I,st . A similar conclusion is obtained 
for the average end-to-end correlation function, which is 
presented in the inset of Fig. 0] 
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FIG. 3: The average string order parameter, O z , the average 
end-to-end correlation function, Ci(L) and the disorder in- 
duced dynamical exponent, 1/z, as a function of the strength 
of disorder. Inset: finite-size dependence of the string order 
parameter and the average end-to-end correlation function 
at the critical point (8 C = 1, open symbols) and in the RS 
phase (S c = 1.5, full symbols) in a log-log plot. The slope 
of the broken lines representing the critical exponents are 
o, 0.39 ± 0.03(0.509]; •, 0.41 ± 0.04[0.382]; 0, 0.69 ± 0.05; 
and ♦, 0.86 ± 0.06(1.0] where in the brackets we presented 
the theoretical RG results; see Table [I] Typical values of the 
error are indicated; otherwise, the error is smaller than the 
size of the symbol. 



O 0.4 




FIG. 4: The average string orderparameter as a function of 
L~ Vst , with rj 3 t — 0.39 as obtained in the inset of Fig. [3] for 
disorder, S = 0.3, 0.5, 0.8, 1.0, and 1.5 from the top to the 
bottom. Inset: the average end-to-end correlation function 
as a function of L~ V1 , 7/1 = 0.69 being the critical decay 
exponent, for the same values of disorder as in the main panel. 
Solid straight lines over the 8 = 1.0 points are guide to the 
eyes. 



C. Critical point and the RS phase 

Our aim with the numerical investigations in this sub- 
section is twofold: first, to check the properties of the 
RS phase, thus to present numerical evidences, and sec- 
ond, and this is numerically more demanding, to try to 
discriminate between the properties in RS phase and at 



the critical point. We start to analyze the finite-size de- 
pendence of the average string order parameter and that 
of the average end-to-end correlation function, which is 
shown in the inset of Fig. [31 at two values of the disorder, 
S = 1 and S = 1.5. The first value should be close to the 
critical point (see Fig. [3J; however, there is certainl;y 
some uncertainty, see the values of Wg in Table ITT1 The 
second value of disorder, S = 1.5, should be deeply in the 
RS phase; however, see the RG estimates in Ref. 0. 

At 6 = 1.5 the decay of the average string order pa- 
rameter, as well as that of the end-to-end correlation 
function is algebraic, and the decay exponents of both 
quantities are in agreement with the theoretical predic- 
tion in the RS phase, as given in Table D For the av- 
erage end-to-end correlation function we have somewhat 
less accuracy, which could be due to similar crossover ef- 
fects as noted for the bulk spin-spin correlation function 
in Ref. The same type of analysis of the results at 
S = 1 give somewhat different results. The decay of the 
average end-to-end correlation function is characterized 
by an exponent, rji — 0.69, which differs considerably 
from the value in the RS phase. On the other hand, the 
decay exponent of the average string order parameter, 
within the error of the calculation, agrees with the value 
found in the RS phase. We note that at the same disor- 
der in the QMC simulation 13 , also the exponent in the RS 
phase is observed. One possible explanation is that 6=1 
is already in the RS phase and therefore we find the cor- 
responding exponent. Anyway, even at the critical point 
one expects strong crossover effects due to the vicinity of 
the RS fixed point, so that probably much larger systems 
are needed to observe the true asymptotic behavior. 

Finally, we compare in Fig. [3] the distribution of the 
gaps at the critical point (a) and in the RS phase (b). 
For both cases the distribution is broadened with in- 
creasing L, which is a clear signal of infinite disorder 
scaling. Indeed, one can obtain a good scaling collapse 
using the form in Eq. I|14|). In the insets we have il- 
lustrated this type of behavior by using the theoretical 
predictions, tp = 1/3 at the critical point and "0 = 1/2 
in the RS phase, respectively. The estimated exponents 
obtained from the optimal scaling collapse are shown in 
Tabled 



IV. CONCLUSION 

The random antiferromagnetic S — 1 chain is a 
paradigm of disorder induced phase transition phenom- 
ena (see also Ref. \34) for which detailed strong disorder 
RG predictions are available. These predictions, how- 
ever, have only been partially verified by numerical calcu- 
lations and even the numerical results are somewhat con- 
flicting. In this paper we have used extensive DMRG cal- 
culations with the aim to clarify the low-energy proper- 
ties of the system with varying strengths of disorder. The 
sizes of the systems we used in the calculation are compa- 
rable with those in previous DMRG studiesjii however, 



InA/L* InA/L 



-6 -3 -8-4 




-25 -15 -5 -35 -25 -15 -5 



InA InA 

FIG. 5: (Color online) Distribution of the gaps in finite sys- 
tems at the critical point, 8 = 1 (a); and in the RS phase, 
S = 1.4 (b). In the insets scaling collapse with Eq. 1141 is 
shown with -0 = 1/3 (a) and ip = 1/2 (b). 

we used a power-law distribution of the couplings in Eq. 
@, which can be more random, than the box distribu- 
tion in Eq. © used previously. We have also calculated 
a quantity, the end-to-end spin-spin correlation function, 
which can be used to locate the borders of the differ- 
ent phases in the system and to obtain an independent 
estimate of the dynamical exponent. Our calculations 
gave further numerical support of the phase diagram pre- 
dicted by the strong disorder RG method and our results 
are basically in agreement with the scenario of disorder 
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induced phase transitions. In the RS phase we made cal- 
culations far from the critical point, which is not possible 
with boxlike distribution of couplings as given in Eq. © 
and obtained estimates for the critical exponents which 
are compatible with the RG predictions. Our results at 
the critical point are less conclusive, which is probably 
due to crossover effects and the inaccurate location of 
this point. For the critical exponent of the end-to-end 
correlation function, 771, and that of the gap scaling, ip, 
numerical estimates at the critical point are clearly dif- 
ferent from that in the RS phase, which are in accordance 
with the RG results. On the other hand, for the average 
string correlation function our numerical results are in 
conflict with the RG prediction. We believe that at this 
point much larger finite systems are necessary to obtain 
a precise numerical estimate and thus to be able to test 
the results of RG predictions. 

We close our paper by mentioning that the present 
day numerical possibilities to explore the properties of 
the random S = 1 chain seem to be exhausted, as far as 
DMRG or QMC methods are considered. Some indepen- 
dent and probably more accurate results can be expected, 
however, by the numerical application of different vari- 
ants of the strong disorder RG method, in particular in 
the vicinity of the critical point and in the RS phase. Re- 
sults obtained in this direction will be published in the 
future^ 
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